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This paper introduces stochastic processes that describe the evo- 
lution of systems of particles in which particles immigrate according 
to a Poisson measure and split according to a self-similar fragmenta- 
tion. Criteria for existence and absence of stationary distributions are 
established and uniqueness is proved. Also, convergence rates to the 
stationary distribution are given. Linear equations which are the de- 
terministic counterparts of fragmentation with immigration processes 
are next considered. As in the stochastic case, existence and unique- 
ness of solutions, as well as existence and uniqueness of stationary 
solutions, are investigated. 

1. Introduction. The aim of this paper is to study random and deter- 
ministic models that describe the evolution of systems of particles in which 
two independent phenomena take place: immigration and fragmentation of 
particles. Particles immigrate and split into smaller particles, which, in turn, 
continue splitting, at rates that depend on their mass. Such a situation oc- 
curs, for example, in grinding lines [1, 23] where macroscopic blocks are 
continuously placed in tumbling ball mills that reduce them to microscopic 
fragments. These microscopic fragments then undergo a chemical process to 
extract the minerals. In such systems, one may expect to attain an equilib- 
rium, as the immigration may compensate for the fragmentation of particles. 
The investigation of existence and uniqueness of such stationary state, as 
well as convergence to the stationary state, is one of the main points of in- 
terest of this paper. It will be undertaken both in random and deterministic 
settings. 

We first introduce continuous-time fragmentation with immigration Markov 
processes. Roughly, their dynamics are described as follows. The immigra- 
tion is coded by a Poisson measure with intensity I[ds) dt, t>0, where / is a 
measure supported on V, the set of decreasing sequences s = {sj,j > 1) that 
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converge to 0. That is, if (s(tj),tj) denotes the atoms of this Poisson mea- 
sure, a group of particles with masses (si(ti), 52(^1), ■ • ■ ) immigrates at time 
ti for each ti > 0. We further impose that / integrates J2j>i{sj ^ !)> which 
means that the total mass of immigrants on a finite time interval is finite 
a.s. The particles fragment independently of the immigration, according to 
a "self-similar fragmentation with index a € M" as introduced by Bertoin 
in [6, 7]. This means that each particle splits independently of others with 
a rate proportional to its mass to the power a and that the resulting par- 
ticles continue splitting with the same rules. Rigorous definitions are given 
in Sections 1.1 and 1.2 below. Some examples of such processes arise from 
classical stochastic processes, as a Brownian motion with positive drift. This 
is detailed in Section 4. 

Let FI denote a fragmentation with immigration process. Our first pur- 
pose is to know whether it is possible to find a stationary distribution for FI . 
Let us mention here that, until now, an equilibrium could only be obtained 
for fragmentation with coagulation processes. See, for example, [4, 14, 16]. 

Under some conditions that depend both on the dynamics of the frag- 
mentation and on the immigration, we construct a random variable Ugtat 
in T> whose distribution is stationary for FI. Let aj be the /-dependent 
parameter defined by 

(1) aj := —supl^a> 0: J s^l{s^>i}/(fis) < oo|. 

When a/ < 0, we obtain that the stationary state Ustat exists as soon as the 
index of self-similarity a is larger than aj and that there is no stationary 
distribution when a is smaller than aj. In this latter case, the particles with 
mass larger than 1, which split slower when a is smaller, do not split fast 
enough to compensate the immigration of large particles, which therefore 
accumulate. Li others words, too many large particles are brought in the 
ball mill which is not able to grind them fast enough. These results are made 
precise in Theorems 7, 8 and 9, Section 2, where we also study whether Ugtat 
is in P, p>0. In addition, the stationary solution is proved unique. 
It is easily checked from the construction of Usjtat that 

Ustat 

as soon as the stationary distribution exists and that this convergence holds 
independently of the initial distribution. One standard problem is to in- 
vestigate the rate of convergence to this stationary state. Our approach is 
based on a coupling method. This provides rates of convergence that differ 
significantly according as a < 0, a = or a > 0: one obtains that the conver- 
gence takes place at a geometric rate when a = 0, at rate when a > 0, 
whereas the rate of convergence depends both on / and a when a < 0. 
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We next turn to deterministic models, namely, fragmentation with im- 
migration equations. Roughly, these equations are obtained by adding an 
immigration term to a family of well-known fragmentation equations with 
mass loss [17, 18, 24]: we consider that particles with mass in the interval 
{x,x + dx) arrive at rate fii{dx) which is defined from / by 

/ fix)Mdx):= / 5:/(s,)/(ds), 
Jo ■'^j>i 

for all positive measurable functions /. Solutions to the fragmentation with 
immigration equation do not always exist. We give conditions for existence 
and then show uniqueness. The obtained solution is closely related to the 
stochastic model {FI{t),t > 0): it is — in a sense to be specified — related to 
the expectations of the random measures X]fc>i ^Fih(t)i * > 0. In this deter- 
ministic setting, one may also expect the existence of stationary solutions. 
Provided the average mass immigrated by unit time is finite, we construct 
explicitly a stationary solution which is proved unique. Note that here the 
hypothesis for existence only involves /, not a, contrary to the stochastic 
case. 

This paper is organized as follows. In the remainder of this section we 
first review the definition and some properties of self-similar fragmentations 
(Section 1.1), then we set down the definition of fragmentation with immi- 
gration processes (Section 1.2). The study of existence and uniqueness of a 
stationary distribution is undertaken in Section 2, where we also give criteria 
for existence of a stationary distribution for more general Markov processes 
with immigration. In Section 3 we investigate the rate of convergence to the 
stationary distribution. Section 4 is devoted to examples of fragmentation 
with immigration processes constructed from Brownian motions with posi- 
tive drift. Section 5 concerns the fragmentation with immigration equation. 



1.1. Self- similar fragmentations. 

State space. We endow the state space 

P = |s = (sj)j>i : si > S2 > • • • > 0, lim Sj = o| 



with the uniform distance 

(i(s, s') := sup |sj — s'A. 

Clearly, as n — > oo, d{s, s*^) — > is equivalent to s" — > Sj for all j > 1, which, 
in turn, is equivalent to Y^j>if{sJ) J2j>ifisj) for all continuous func- 
tions / with compact support in (0, oo). Hence, P identifies with the set of 
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Radon counting measures on (0,oo) with bounded support endowed with 
the topology of vague convergence through the homeomorphism 

With a shght abuse of notation, we also call s the measure X]j>i'5sjl{sj>o}- 
It is then natural to denote by "s + s'" the decreasing rearrangement of the 
concatenation of sequences s,s' and by (s, /) the sum X]j>i /(si)l{sj>o}- 
More generally, we denote by "X]i>iS*" the measure J2i>iJ2j>i^s'^.^{s^.>o}- 
This point measure does not necessarily correspond to a sequence in V, but 
when it does, it represents the decreasing rearrangement of the concatenation 
of sequences s^, s^, . . . . 

For all p>0, let P be the subset of D of sequences si > S2 ^ ■ ■ ■ ^ such 
that J2j>i'^^ < oo, endowed with the topology of V. When p = 0, we use 
the convention 0*^ = 0, which means that is the space of sequences with 
at most a finite number of nonzero terms. Let also T>i be the subset of T> of 
sequences such that J2j>i ^ 1- Clearly, P C F when p <p' and Vi Cl^. 
At last, set := (0, 0, . . ."). 

Self-similar fragmentations. 

Definition 1. A standard self-similar fragmentation {F{t),t > 0) with 
index a G M is a Pi-valued Markov process continuous in probability such 
that: 

. F(0) = (1,0,...), 

• for each to > 0, conditionally on -F(io) = (■^i, S2, • • • )' the process + 
to),t > 0) has the same law as the process obtained for each t > by 
ranking in the decreasing order the components of sequences siF^^^sft), 
S2F^'^\s^t), where the F(j)'s are independent copies of F. 

This means that the particles present at a time to evolve independently 
and that the evolution process of a particle with mass m has the same dis- 
tribution as m times the process starting from a particle with mass 1, up 
to the time change t i-^ tm°^ . According to [3] and [7] , a self-similar frag- 
mentation is Feller — hence, possesses a cadlag version which we shall always 
consider — and its distribution is characterized by a 3-tuple (a,c, z^): a is 
the index of self-similarity, c > an erosion coefficient and v a dislocation 
measure, which is a sigma-finite nonnegative measure on T>i that does not 
charge (1,0,...) and satisfies 



f (1 — si)v{ds) < oo. 
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Roughly speaking, the erosion is a deterministic continuous phenomenon 
and the dislocation measure describes the rates of sudden dislocations: a 
fragment with mass m splits into fragments with masses ms,s € Pi, at rate 
m"z/((is). In case z^(Pi) < oo and c = 0, this means that a particle with mass 
m splits after a time T with an exponential law with parameter m'^viVi) 
into particles with masses ms, where s is distributed according to z^(-) /uiVi) 
and is independent of T. For more details on these fundamental properties 
of self-similar fragmentations, we refer to [3, 6, 7]. 

Definition 2. For any random u G D, a fragmentation process (a, c, i^), 
starting from u, is defined by 

(2) F'^-\t):=Y.{u,F^'\u'-t)), t>0, 

i>i 

where the F^^^'s are i.i.d. copies of a standard (a, c, z^)-fragmentation F, 
independent of u. 

Clearly, F^^\t) G V for alH > and, according to the branching property 
of F, F^") is Markov. It is plain that such a fragmentation process converges 
a.s. to as t ^ oo, provided z^(I'i) 7^ 0. We shall denote in the sequel by 

F^^\t) > • • • > -F^"^(t) > • • • the components of the sequence F^^'^{t). 

We now review some facts about standard (a, c, z^)-fragmentations that 
we will need. In the remainder of this section, F denotes a standard (a, c, v)- 
fragmentation. 

Tagged particle. We are interested in the evolution process A of the mass 
of a particle tagged at random in the fragmentation. To construct this pro- 
cess, we recall that one may always suppose that F is built from some family 
{G{t),t > 0) of nested open sets of (0, 1) so that F{t) is the ordered sequence 
of lengths of the interval components of G(t), t > (see [3, 7]). Let then U 
be uniformly distributed on (0, 1), independent of G, and call \{t) the length 
of the interval component of G{t) containing U. When such interval does 
not exist, set \{t) := 0. The main point of interest of such approach is that 
the distribution of A is well known. 

First, when a = 0, Bertoin [6] shows that A =^ exp(— ^(•)), where ^ is a 
subordinator (i.e., a right-continuous increasing process with values in [0, co] 
and with stationary and independent increments on the interval {t:^{t) < 
00}), with Laplace exponent 4> given by 

(3) ^{q) := c{q + l)+ [ f 1 - ^ s]+A u{ds), q > 0. 

"'^i V j>i ) 

We recall that (\) characterizes ^, since i?[exp(— (/^(t))] = exp(— t(^(g)) for all 
t,q>0 (for background on subordinators, we refer to [5], Chapter III). When 
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c > or i^{J2j>i < 1) > Oi one sees that the subordinator ^ is killed at rate 
k = (/)(0) > 0: that is, there exists a subordinator with Laplace exponent 
(f) = (p — k and an exponential r.v. e(k) with parameter k, independent of ^, 
such that 

Cit) =?(*)l{t<e(fe)} + 00l{t>e(fe)} 

for ah t > 0. 

When a E M, Bertoin [7] shows that A exp(— ^(p(-))), where ^ is the 
same subordinator as above and p is the time-change 

(4) :=inf|ii>0:^ exp(a^(r))dr >t|, t > 0. 

On the other hand, the construction of A implies that, conditionally on F, 
X{t) = Fk{t) with probability Ff^{t), k> 1, and that A(t) =0 with probability 
1 - J2k>iPk{t)- Hence, 

(5) 5:/(F,(t)) = i^[/(exp(-^(p(t))))exp(e(Ht)))|i^] 

k>l 

for every positive measurable function / supported on a compact subset of 
(0,oo) (with the convention x cxo = 0), and, in particular. 



(6) E 



k>l 



E[fiexpi-apm)expiapm]- 



Formation of dust when a < 0. When the index of self-similarity a is 
negative, for all dislocation measures u, the total mass J2k>i ^kit) of the 
fragmentation F decreases as time passes to reach in finiFe time even if 
there is no erosion (c = 0) and no mass is lost within sudden dislocations 
{i^Cl2j>iSj < 1) = 0). This is due to an intensive fragmentation of small 
particles which reduces macroscopic particles to an infinite number of zero- 
mass particles or dust. To say this precisely, introduce 



(7) 



r:=inf|i>0:5^Ffc(t)=ol 
I k>l J 



the first time at which the total mass reaches 0. According to Proposition 14 
in [18], there exist C,C' some positive finite constants such that, for any 
t>0, 

(8) P{T>t)<Cexp{-C't^), 
where F is a (c, i^)-dependent parameter defined by 

(1 — A)~^, when (j){q) — cq 

(9) r := varies regularly with index 0<A<lasg^oo, 
1, otherwise. 
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Note that E[t] < oo. This phenomenon of formation of dust does not occur 
when a > 0: if no mass is lost by erosion or within sudden dislocations, then 
Y.k>i Pkit) = 1 a.s. for all t > 0. 

1.2. Fragmentation with immigration processes. As said previously, the 
immigration and fragmentation phenomena occur independently. The im- 
migration is coded by a Poisson measure on x [0, oo) with an intensity 
I{ds) dt such that 




and we call such measure / an immigration measure. The hypothesis (HI) 
implies that the total mass of particles that have immigrated during a time 
t is almost surely finite (for an introduction to Poisson measures, we refer 
to [21]). On the other hand, the particles fragment according to a self-similar 
fragmentation (a,c, z^). 

Definition 3. Let u be a random sequence of T) and let ({s{ti),ti),i > 
1) be the atoms of a Poisson measure with intensity I{ds)dt independent 
of u. Then, conditionally on u and ((s(ti), tj), i > 1), let > 1, 

be independent fragmentation processes {a, c, i^) starting, respectively, from 
u, s(ti), 3(^2), • • • • With probability one, the sum 

(t) := (t) + (t - ti) 

belongs to P for all t > 0, and the process F/^"^ is called a fragmentation 
with immigration process with parameters {a,c,u,I) starting from u. 

One may be troubled by conditioning on the value of {{s{ti),ti),i > 1), as 
it may have probability. If so, note that the family > 1, is actually 

constructed from the Poisson measure ((s(tj), tj), i > 1) and an independent 
family F^^'^\i, j > 1, of i.i.d. standard (q, c, i^) -fragmentations, through the 
formula F«'^^\t) =j:ij>^Sj{t,)F(''^\{.s,{ti)rt). 

The reason why Et,<t ^^^^^^''^^-^i) 6 ^ a.s. is that Et,<t Ej>i Sj{ti) < 00 
[by hypothesis (HI)] and then that Et,<t ^^'^**"(i-*i) G ^^ since Efc>i ^f^*'^^(*- 
ti) < Ei>i Sj{ti). Note also that when p > 1, G F as soon as u€ 1^. 

In this definition, the sequence u represents the masses of particles present 
at time and at each time tj > 0, some particles of masses s{ti) immigrate. 
At time t, two families of particles are then present: those resulting from the 
fragmentation of u during a time t and those resulting from the fragmenta- 
tion of s(tj) during a time t — ti, ti<t. 
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It is easy to see that the process FI^^^ is Markov with the Feller property 
(cf. the proof of Proposition 1.1 in [3]). Hence, we may and will always 
consider cadlag versions of F7("). 

In the rest of this paper, we denote by FI a fragmentation with immigra- 
tion {a, 0,1^,1) (without any specified starting point) and we always exclude 
the trivial cases u = or 1 = 0. 

Remark. One may wonder why we do not more generally consider some 
fragmentation with immigration processes with values in TZ, the set of Radon 
point measures on (0, oo). Indeed, for all (random) u G 7^ and all t > 0, it is 
always possible to define the point measure 

(10) F/(")(t):=F(")(t) + ^F(^(*'»(t-ti), t>0, 

where F^'^\t) is defined similarly to (2) and is independent of f(®(*')),z > 1, 
some independent fragmentations (a,c, z^) starting, respectively, from s(ti), 

s{t2), The sum involving the terms F(®(*^))(t — ti), ti < t, is in T>, as 

noticed in the Definition 3 above. The issue is that, in general, starting from 
some u € TZyD, the measures F^^\t) do not necessarily belong to TZ, as 
the masses of the initial particles may accumulate in some bounded interval 
(a, b) after fragmentation. 

As an example, starts from u = J2i>i ^-^d fix a > 0. For each i, tag a 
particle at random in the fragmentation issued from the particle with mass 
I, as explained in the previous section. At time i, this tagged particle is 
distributed as i exp(— ^(*)(/9(*)(i"i))), where the .^(*^'s are i.i.d. subordina- 
tors with Laplace exponent (3) and p^'^^ the corresponding time changes (4). 
According to the Borel-Cantelli lemma, the number of tagged particles be- 
longing to some interval (a, h) at time t is then a.s. infinite [and, there- 
fore, F(")(t) ^7^] as soon as ^i>iP{a < i exp(— ^(p(i"t))) < 6) = oo. In [9], 
Bertoin and Caballero show that for most of subordinators (and, therefore, 
for most of dislocation measures) i exp(— ^(/3(i"t))) has a nontrivial limiting 
distribution as i — > oo when a > 0. In such cases, the above sum of proba- 
bilities is infinite for some well-chosen intervals (a, h) and then F^") {t) ^ TZ. 

That is why we study fragmentation with immigration processes on P. 
However, in Section 5, we shall use some of these measures 
and we give (Proposition 15) some sufficient conditions on u and a for 
-F(")(t) [equivalently, FI^^\t)] to be a.s. Radon at fixed time t. These con- 
ditions do not ensure that the process F/(") is 7^-valued, as we do not know 
if a.s., for all t, FI^'^\t) £TZ. 

2. Existence and uniqueness of the stationary distribution. This section 
is devoted to the existence and uniqueness of a stationary distribution for 
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FI and to properties of the stationary state, when it exists. We begin by 
estabhshing some criteria for existence and uniqueness of a stationary distri- 
bution, which are available for a class of Markov processes with immigration 
including fragmentation with immigration processes. This is undertaken in 
Section 2.1 where we more specifically obtain an explicit construction of a 
stationary state. We then apply these results to fragmentation with immi- 
gration processes (Section 2.2). 

From now on, for any r.v. X, C{X) denotes the distribution of X. 

2.1. The candidate for a stationary distribution for Markov processes with 
immigration. Recall that IZ denotes the set of Radon point measures on 
(0, oo) and equip it with the topology of vague convergence. We first con- 
sider 7^- valued Markov processes with some superposition property and then 
extend the results to a larger class of Markov processes. 

Let X be an 7?.-valued Markov process that satisfies the following superpo- 
sition property: for all u, v G 7^, the sum of two independent processes X^") 
and X^^^ starting, respectively, from u and v is distributed as X^^~^'^\ A 

moment of thought shows that this is equivalent to X]i>i X^^''^ ^= X^^'^i " 
for all sequences (u*, i > 1) such that Y.i>i £ 7^ a.s., where , X'^'^'\ . . . 

are independent processes, starting, respectively, from u^,u^, Consider 

then J, a nonnegative cr-finite measure on 7^, and let {{s{ti),ti),i > 1) be 
the atoms of a Poisson measure with intensity I{ds) dt,t > 0. Conditionally 
on this Poisson measure, let be independent versions of X, start- 

ing, respectively, from s(ti), s(t2), • . • . In order to define an X-process with 
immigration, we need and will suppose in this section that a.s. 

J2x^<^^^\t-ti)en forant>0. 

ti<t 

In particular, this holds when X is a fragmentation process and I an im- 
migration measure, as explained just after Definition 3. More generally, still 
supposing that X is a fragmentation, one easily checked that it holds as 
soon as / integrates Ijs^^g} for all e > 0, that is, as soon as the number of 
particles of mass larger than e immigrating in finite time is finite. 

Definition 4. For every random u G 7^, let X^^^ be a version of X start- 
ing from u and consider {{X^^^'"^^\vi),i > 1) a version of ((X(^(*»)),ti),i> 1) 
independent of X^^\ Then, the process defined by 

(11) X/(")(t):=X(")(t) + ^X(''(''>))(t-^;i), t>0, 

Vi<t 

is an 7^-valued Markov process and is called X-process with immigration 
starting from u. 
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We point out that the Markov property of XI results both from the 
Markov property and from the superposition property oi X . A moment of 
reflection shows that the law of the point measure 

(12) Ustat:=E^^'^*'"(*^) 

ti>0 

is a natural candidate for a stationary distribution for XI [in some sense, 
it is the limit as t ^ oo of XI^^\t)], provided that it belongs to TZ. The 
problem is that it does not necessarily belong to TZ, as the components of 
Ustat may accumulate in some bounded interval {a,b). 

Lemma 5. (i) //Ustat £ TZ a.s., then the distribution £(Ustat) is a sta- 

f \ p 

tionary distribution for XI and for any random u£TZ such that X^^> (t) — > 
as t ^ oo, 

Ustat ast^oo. 

(ii) If -P(Ugtat ^TZ) > 0, then there exists no stationary distribution for 
XI and if P(Ustat ^ ^) > 0, then there exists no stationary distribution on 
V for XI. 

Proof, (i) Assume Ustat £ IZ a.s. and consider a version 

X/(Ustat) of 

the X-process with immigration starting from Ustat- We want to prove that 
X/(UBtat)(i) 'if u^t^t for every t > 0. So fix t > 0. By definition of XI and 
using the Markov and superposition properties of X, we see that there exists 
{{X'^<''-^\vi),i> 1) an independent copy of > 1) such that 

ti>0 Vi<t 

By independence of ((r(f j), Uj), i > 1) and {{s{ti),ti),i > 1), the concatena- 
tion of 

{{r{vi),t-Vi),Vi<t) and {{s{ti),ti + t),i>l) 
has the same law as ((s(ij), tj), i > 1). Hence, 

^j(U..O(i)l|?'^X(^(*.))(i^)= Ustat. 
ti>0 

Similarly, one obtains that, for all t>0, 

(13) X/(")(i)^= X(")(t) + E^^''^^'^^('^i), 

Vi<t 
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where {{X^''^''''>\vi),i > 1) is distributed as {{X^<^''^\ti),i > 1) and is inde- 
pendent of X^^\ Suppose now that (t) as t ^ oo. Clearly, 

Vi<t fi>0 

and, therefore, 

(t) + ^^'^"'^^ (^i) ^ E ^^'^'"^^ (^i) as t ^ oo. 
Since the limit here is distributed as Ustat and since (13) holds, one has 

X/(")(i)'^Ustat. 

(ii) Suppose that there exists a stationary distribution £stat- Our aim is 
to show that P(Ustat ^ 7^) = 0. To do so, let Xl'^^-^^^^ be an X -process with 
immigration starting from an initial sequence distributed according to £stat • 
Replacing u by X/(^=*^*)(0) in (13), we get 

J^j(£stat ) (0) ^(X/(^Bta.) (0)) ^ ^ ^(s(tO) 

Introduce then, for any < a < 6 < oo, the event 

i?,,,:=|^(X(^(*'»(tO,l{a,fe))=oo| 
Ui>0 J 

and fix some > 0. The identity in law obtained above yields 
P((X/(^--)(0),l(,,,))<iV) 

<p('^(X(^(*>))(ti),l(,,,))<iv') 
\ti<t I 

\u<t J 

The first probability in this latter sum converges to as t — > oo by definition 
of Ea^b and, therefore, 

P((X/(^stat)(o), < AT) < P{n\Ea,b) ViV > 0. 

Letting — > oo, we get P{0,\Eafi) = 1 (because £stat is supported on TZ) 
and then P{Ea,b) = 0. This imphes that P(Ustat ^ 7^) = 0. 

Now, replacing 7^ by P and Ea^b by Ea^oo, we obtain similarly that 
-P(Ustat ^V) =0 as soon as there exists a stationary distribution Cgtut such 
that £stat(^^) = l. □ 
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Let us now extend these results to Markov processes that take values in 
some a-compact space E and that do not necessarily satisfy the superpo- 
sition property. In order to introduce some immigration and some super- 
position property, we will work on SJl^;, the set of point measures on E: if 
m € TXe, either m = J2i>i^x(^^ some sequence {x^^\i > 1) of points of 
E, or m = 0, where o is the trivial measure: o{E) = 0. The subset of mea- 
sures of QJt^; that are Radon is denoted by 971^^^°'^ and is equipped with 
the topology of vague convergence. Consider then /, a nonnegative a-finite 
measure on E, and {X(t),t > 0), a Markov process with values in E. For 
any m = J2i>i "^iCoG ^E, set 

^(-)(t):=5:<5^(,„,^^^, t>0, 

i>l 

where . . . are independent versions of X, starting, respec- 

tively, from , 2;(2) , . . . . If m = 0, (t) := 0, V t > 0. 

We now construct some X-process with immigration. Let ro be a random 
element of dJl^^'^°'^ and {{x{ti),ti),i > 1) be the atoms of a Poisson measure 
with intensity I{ds)dt, t>0, independent of m. Conditionally on this Pois- 
son measure and on m, let X^^^ and > 1, be independent versions 
of X starting, respectively, from m, Sx(ti)i^x{t2)j ■ ■ ■ ■ Define then 

:= + J2 ^^^'''■'^^\t -ti), t > 0, 

ti<t 

and suppose that a.s., for all t > 0, AfJ^"') G 9Jt|^'^°'^. Then is Marko- 

vian and called A'-process with immigration starting from m. Introduce next 
the point measure 

ti>0 i>l 

With the same kind of arguments as above, one obtains the following result. 

Lemma 6. (i) AssumeUstat G a.s. Then the distribution C{Usta.t) 

is a stationary distribution for X"] and X'i'^'^\t) Z^stat soon as X'^'"^\t) ^ 
as i —> oo. 

(ii) // P(Usta.t ^ OJt^;^'^"'^) > 0, there exists no stationary distribution for 
X3. 

2.2. Conditions for existence and properties of FIs stationary distribu- 
tion. Up to now, / is an immigration measure as defined in Section 1.2, 
that is, / satisfies hypothesis (HI). Let FI denote a fragmentation with im- 
migration (a, c, 1^,1). By definition, the fragmentation process satisfies the 
superposition property and, for every uGD, F^'^\t) as t — > oo. Then 
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the results of Lemma 5 can be rephrased as follows: if ((s(tj), tj), i > 1) are 
the atoms of a Poisson measure with intensity I{ds) dt and if condition- 
ally on this Poisson measure, i?'(®(*i))^^{s(i2))^ _ _ _ ^j^g independent (a,c,z^)- 
fragmentations starting, respectively, from s(ti), s(t2), • • • , then there is a 
stationary distribution for the fragmentation with immigration {a,c,u,I) if 
and only if 

Ustat = E^^'^*^"(*OeP a.s. 

In this case, 

i7/(u)(i)i-'Ustat ast^oo 

for all u G P and, therefore, vC(Ustat) is the unique stationary distribution 
for FI. The point is then to see when Ugtat belongs to P and when it 
does not. The results are given in Section 2.2.1 where we further investigate 
whether Ustat is in F or not, p>0. This is particularly interesting when 
Ustat £ a-s.: then the total mass of the system converges to an equilibrium, 
which means that the immigration compensates the mass lost by formation 
of dust (when a < 0), by erosion or within sudden dislocations. When Ugtat G 
V a.s., we also investigate the behavior of its small components. The proofs 
are detailed in Section 2.2.2. 

2.2.1. Statement of results. Let F denote a standard (a, c, z^)-fragmentation. 
In the statements below, we shall sometimes suppose that 

(H2) c = 0, ui^Sj <l\=0 and / ^\ln{sj)\sji^{ds) < oo 
\j>i J -^^1 j>i 

or 

(H3) $0 < r < 1 : Fi{t) £ {r"- ,n gN} Vt>0,i>l, and (H2) holds. 

In terms of ^, the subordinator driving a tagged fragment of F, the hypoth- 
esis (H2) means that £'[^(1)] < oo. We shall also use the convention F = P 
when p<0. 

We now state our results on the existence of a stationary distribution; 
they depend heavily on the value of the index a. 

Theorem 7. Suppose a < 0. 

(i) If either III J2j>iSj'^'^{sj>i}I{ds) <^ or Jii s{'^ In sil{s^>iy I (ds) < 
oo, then the stationary state Ugtat £ o,-S- for all p> 1 + a. 

(ii) There exists no stationary distribution when J^i sj~"l{sj>i}.I(ds) = oo. 



Theorem 8. Suppose a = 0. 
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(i) // Jii lnsil|s^>i}/((is) < cxD, then, with probability one, Ugtat £ for 
all p> 1 and does not belong to when c = and i^(X]j>i '^j < 1) = 0. 

(ii) There exists no stationary distribution when /^i lnsil{^^>i|/((is) = oo 
and (H2) holds. 

Theorem 9. Suppose a > 0. // /^i sf l{5^>i}/((is) < oo for some e > 0, 
then Ustat £ 0,-S- for p large enough and if (H3) holds, then Ustat ^ 
a.s. More precisely, for every 7 > 0; 

(i) /ii Ej>i'5jl{s3>i}^(c^s) < 00, i/ien Ustat G /or all p > 1 + 
a/(7Al), 

(ii) s?l{,,>i|/(ds) =00 and (H3) /loWs, i/ien Ustat ^ a.s. 

When —1 < a < 0, the result of Theorem 7(i) can be completed (see the 
remark following Proposition 10 below): in most cases, either Ustat ^ 1^'^°' 
a.s. or both events {Ustat = 0} and {Ustat ^ 1^'^°'} have positive probabili- 
ties. 

It is interesting to notice that the above conditions for existence or absence 
of a stationary distribution depend only on a and /, provided hypothesis 
(H3) holds. Indeed, recall the definition (1) of aj and let then a vary. Ac- 
cording to the above theorems, the values a = aj and a = — 1 are critical. 
Provided aj <0, the stationary distribution exists when a > aj and does 
not exist when a < a/. Moreover, the stationary state Ustat is a.s. com- 
posed by a finite number of particles as soon as aj < a < —1, whereas when 
a > — 1, Ustat ^ 1^^°' with a positive probability (which equals 1 when a > 
and depends on further hypothesis on / and a when —1 < a < 0, see the 
forthcoming Proposition 10). 

Let us try to explain these results. By the scaling property of fragmen- 
tation processes, particles with mass > 1 split faster when a is larger. This 
explains that, when a is too small, some particles may accumulate in inter- 
vals of type (a, 00), a > 0, which implies that Ustat ^ T^- For a large enough, 
particles with mass > 1 become rapidly smaller, but particles with mass < 1 
split more slowly when a is larger. Therefore, small particles accumulate and 
Ustat ^ l-^ when p is too small. Moreover, the smallest p such that Ustat G l-^ 
increases as a increases. When a < —1, it is known that small particles are 
very quickly reduced to dust (see, e.g.. Proposition 2 in [8]). This implies 
that Ustat G J-^, provided it belongs to V. 

Small particles behavior. Suppose that — 1 < a < and J^i J2j>i ^J°''^{sj>i} ^ 
I{ds) < 00, so that Ustat S V a.s., according to Theorem 7(i). Consider then 
the random function 

ei-^-Ustat(e) :=Ustat([e,oo)), 
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which counts the number of components of Ustat larger than e. We want to 
investigate the hmiting behavior of UgtatC^) as e — > 0. In that aim, we make 
the following technical hypothesis: 

/ ^ s\^°'sjv{ds) < oo and 
i>i>i 

(H4) 

/ (1 — si)^z/((is) < oo for some 6* < 1, 

as well as hypothesis (H3). Note that the first integral involved in (H4) is 
finite as soon as a > — 1 and z/(sAr > 0) = for some integer N >2, because 
then /^^ E,>.>i sl+'^s.uids) < (N - 1) J^^ (1 - s.Hds). 

Proposition 10. Under the previous hypotheses: 

(i) if JiiJ2j>i ■^J°''^{sj<i}H'^^) < there exists a finite r.v. X, < 
P{X = 0) < 1, inch that 

Ustat (e)e'+" a.s., 

(ii) if Jii s5""ljs^,<i}/((is) = oo, one has liminfe^o e^^°Ustat(e) > a.s. 

In particular, this implies that P(Ustat ^ l^^"') = 1 when the assumption 
of the second statement is satisfied. This is not true when the assumption of 
the first statement holds: in such case, < P(Ustat = 0) < P(Ustat S /-^+") < 
1 [see the proof of (i) for the first inequality]. 

When a > or a < —1, some information on the behavior of Ustat (s) as 
£ — >■ can be deduced from Theorems 7, 8 and 9. Thus, Ustat (0) < oo a.s. as 
soon as aj < a < —1. To obtain some information when a > 0, first notice 
that, when Ustat G ^, /(o,oo) ^^'^stAdx) < oo 44> /(q -^^ UstatCa^^^^) dx < oo, by 
integration by parts. Combined with Theorem 9, this implies, when a > 0, 
that if /;i J2j>i sjl{s^.>i}/((is) < oo, then liminfe^oe^Ustat(e) = for all p > 
l + a/(7Al), whereas if /^i s]'l|^^>;^|/((is) =oo, limsup^^Q ^^Ustat(e) =oo 
for all p < 1 + a/(7 A 1) . The behavior near of Ustat (e) is then strongly con- 
nected to the immigration /. Similarly, when a = and when there is a sta- 
tionary distribution, one deduces from Theorem 8 that lim inf e_>o e^Ustat (s) = 
for all p > 1, and that lim sup^^Q e^Ustat (s) = oo for all p < 1, provided 
c = i^(Ei>iSi < 1) = 0. 

Remark. It is possible to show that Ustat G a.s. as soon as /^i J2j>i ^ 
l{s^>i}/(ds) < oo and that P(Ustat ^ '7^) > as soon as q > —1, /^i sj~°l{sj>i} x 
I{ds) = oo and hypotheses (113) and (114) hold. The first claim can be proved 
by using some arguments of the proof of the forthcoming Proposition 16 and 
the second claim is a consequence of Theorems 4(i) and 7 of [19], which are 
also used below to prove Proposition 10. 
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2.2.2. Proofs. Let F be a standard (a, c, i^)-fragmentation and for every 
p G M and t>0, define 

M(p,t):=5](Ffc(t)fl|^^(,)>0}, 

fc>i 

which is a.s. finite at least when p>l (since it is bounded from above by 1). 
That Ustat belongs to some P-space is closely related to the behavior of the 
function 1 1-^ M{p,t). Indeed, 

Ustat = EE 

where the F^^'^^^s, i,j>l, are i.i.d. copies of F, independent of {{s{ti),ti),i > 
1). Then U^tat £P ^M{p)< oo with 

M(p) = / xPVsMdx) 

i(0,cx>) 

= EE«'(*^)^^'''^(^''^i(*^)*^)iK(t.)>o}> 

where the M(*'-')(p, •)'s, i,j > 1, are i.i.d. copies of M{p,-), independent of 
{{s{ti),ti),i > 1). Using the tagged particle approach as explained in Sec- 
tion 1.1, one obtains the following results on M{p, •). 

Lemma 11. (i) Suppose a<0. Then ex.p{Xt)E[M{p,t)]dt < oo as 
soon as p>l + a and A < (j){p — 1 — a). In particular, E[M{p,t)] < oo for 
a.e. t>0 as soon as p>l + a. 

(ii) Suppose a > 0. Then for every rj > and every p>l, there exists a 
random variable F)(r^^p) with positive moments of all orders such that 

M{p, t) < L>(^^p)t-(P~i)/("+'') a.s. for every t > 0. 
Consequently, E[M {p, t)] dt < oo when p> 1 + a. 

Bertoin (Corollary 3 in [8]) shows that when a > and p> 1, the pro- 
cess tiP~^)/<^M{p,t) converges in probability to some deterministic limit as 
t — > oo, provided the fragmentation satisfies hypothesis (H3). See also Bren- 
nan and Durrett [11, 12], who prove the almost sure convergence for binary 
fragmentations (z/(si -|- S2 < 1) = 0) with a finite dislocation measure. 

Proof. We use the notation introduced in Section 1.1. 
(i) According to (6), 



E[M{p,t)]=E[exp{{l-p)apm'i-{t<D}], 
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where D = inf{t : p{t) > e{k)}. Therefore, 
exp{Xt)E[M{p,t)]dt 

r 

(k) 

exp(Ap"^ (t)) exp((l -p + a)C(t)) dt 
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(14) 



E 



E 



exp(Ai)exp((l-:pK(p(t)))(it 

e(fc) 



using for the last equahty the change of variables 1 1— > p[t) and that, by 
definition of p, exp{a(,{p{t))) dp{t) = dt on [0,D). The function p~^ denotes 
the right inverse of p and, clearly, p~^{t) < t since a < 0. When p>l + a, 
this leads to 

exp{Xt)E[M{p,t)]dt 



< 



exp{—4>{p — 1 — a)t) dt 



e(fc) 



exp((A — — 1 — a))t) dt 



if A < 0, 
if A > 0, 



and in both cases, the integral is finite as soon as A < (l){p — 1 — a) = (p{p 
l-a) + k. 

(ii) Fix a > 0, p >1 and rj > and recall that, according to (5), 

M{p,t) = E[eM-{p-mp{tm{t<D}\F]. 

Since is increasing, one has 

rp{t) /-oo 
p{t) exp{—r]^{p{t))) < / exp{—ri^{r)) dr < / exp{—r]^{r))dr:=D 
Jo Jo 

And, on the other hand, for t < D, 

rp(t) 

exp(a^(r)) dr < p{t) exp(a^(p(t))). 



Combining these inequalities, we obtain exp(— (a + i])C{p{t))) < t ^^(,;) for 
ah t < D. Hence, M{p, t) < t-(P-i)/(°+^)L'(^,p) , where D^^^p) := 
Carmona, Petit and Yor [13] have shown that Z)(^) has moments of all pos- 
itive orders, which, by Holder's inequality, is also true for D^^py □ 

We now turn to the proofs of Theorems 7, 8 and 9. 

Proof of Theorem 7. (i) Fix p>l + a and split M{p) into two sub- 
sums: 
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and Msup(p) = M{p) — Minf(p). One has 

E[MM]=l^ (^^^^spi|,^,<i}^/(ds) X E[M{p,t)]dt 

and both of these integrals are finite according to hypothesis (HI) and 
Lemma 11 since p > 1 + a. It remains to show that Msup(p) < oo when 
/ integrates Ej>i s7°l{s,>i} or s^"" lnsil{s^>i}. 

Suppose first that /^i X]j>i s~j°''^{sj>i}I{d^) < and let r^*'-'^ be the first 
time at which the fragmentation is entirely reduced to dust. Equiv- 

alently, the first time at which M^^'^^ reaches 0. If the number of 

pairs such that s^{ti)ti < r^*'-'^ and Sj(ti) > 1 is finite, then the sum 
Msup(p) is finite because it involves at most a finite number of nonzero 
M^''^\p,s'^{ti)ti) [which are a.s. all finite according to Lemma ll(i)]. To 
prove that this is the case, we use the theory of Poisson measures. Since the 
r.v. T^'^'^\ i,j > 1, are i.i.d., the measure 

i>l ^ 

is a Poisson measure with intensity m defined for any positive measurable 
function / by 



/•oo POO r r / 

/ f{x)m{dx)= / E fit-^ sup (r(i'^)s7 

Jo Jo Jl^ I V j:s,>l ■' 



I{ds) dt. 



The integral m{dx) is bounded from above by E[t*'^'^^] /^i J2j>i ^J°''^{sj>i} ^ 
I{ds), which is finite by assumption on / and since E'fr^^'^^] < oo [by (8)]. 
This implies that a.s. there are only a finite number of integers i > 1 such 
that sup j . g^^t^^-^i{T^^'^^ sj'^ (ti)) > 1. For each of these i, there is at most 
a finite number of integers j > 1 such that Sj{ti) > 1. Hence, the number of 
pairs such that s^{ti)ti < r^*'-'^ and Sj{ti) > 1 is indeed a.s. finite. 

Assume now that /^i lnsil{si>i}-^(f^s) < oo. For any a > 0, the number 
of integers i > 1 such that ati < s^°'{ti) ln(si(ij)) and si(tj) > 1 is then a.s. 
finite. The sum Ms^p{p) is therefore finite if 

is finite for some (and then all) a > 0. The expectation of this latter sum is 
bounded from above by 

^ /- El{«.>i}^(^^) / exp{at{p-a))E[M{p,t)]dt, 



EQUILIBRIUM FOR FRAGMENTATION WITH IMMIGRATION 19 

which is finite for a sufficiently small, according to Lemma ll(i). Hence, 
Miup(?') < oo a.s. 

(ii) Suppose /^i sj~"l{s-^>i}/((is) = oo and let t^'^^ := mi{t > : Fl^'^\t) < 
1/2} be the first time at which all components of are smaller than 1/2, 
i>l. Note that S[t-[^2^] > since F^^'^^ is cadlag. The measure 

i>l:si(t,)>l ^'^ 

is a Poisson measure with intensity m' given by 

f{x)m'{dx) = j^^E[f{s^^r\^lf^^)]l{,^^^^I{ds)dt. 

By assumption on / and since ii'fr^yg'^^] > 0, the integral /^°° m'{dx) is infinite 

and, consequently, the number of integers i such that t^'^^ > Si{ti)ti and 

si{ti) > 1 is a.s. infinite. For those i, si{ti)F^'^\si{ti)ti) > 1/2 and, there- 
fore, Ustat contains a sequence of terms all larger than 1/2, which implies 
that it is not in T> a.s. □ 



Proof of Theorem 8. (i) The second part of the proof of Theorem 7(i) 
(replacing there a by 0) shows that Ustat G rip>i when /^i ln(si)l|<,j>i}/((is) < 
oo. Now, if c = and i^{J2k>i Sfc < 1) = 0, the sum M(l) equals J2i>i J2j>i 
which is clearly a.s. infinite since 1^0. 

(ii) Assume that J^i ln(si)l{5j>i}/((is) = oo and £^[^(1)] < oo. For each 

i > 1, let exp(-^(*'i)(-)) denote the process of masses of the tagged particle 
in the fragmentation F^*'^) . To prove that Ustat ^ T^, it suffices to show that 
its subsequence {si(ti)exp(— ^(*'-^)(ij)),i > 1}-^ ^ V. The components of this 
sequence are the atoms of a Poisson measure with intensity m" given by 

/•oo rOG r 

/ f{x)m"{dx)= E[f{siexpi-m)]Iids)dt. 
Jo Jo Jl^ 

Take then a > i?[^(l)]. Since S,{t)/t —> £^[^(1)] as t — > oo, there exists some 
to such that P(^(t) < at) > 1/2 for t > to- Then 

/ m"{dx)= / / P{i{t)<\usi)I{ds)dt 
Ji Jo Ji^ 

> / / P{^{t)<at)dtl{ds) 

Jl^ Jto 

- ^^^(a"^lnsi - to)l{a-ilnsi>to}^('^s) 
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and this last integral is infinite by assumption. Hence, J2i>i ^si(i,)exp(-5{».i)(t,)) ^ 
D a.s. and a fortiori Ustat ^ ^ a-S- □ 

Proof of Theorem 9. Fix p>l + a. According to the Campbeh 
formula for Poisson measures (see [21]), the sum M{p) is finite if and only if 



(15) 



Ji 



1-exp -^s^M(1'J)(p, sp) 



I{ds) dt < oo. 



(i) We first prove assertion (i) and that Ugtat G a-S- for p large enough 
when / integrates sf l|^^>i}. Suppose p> 1 + a and note that the integral 
(15) is bounded from above by 

//•oo 
Vspi{,^<i}/(ds) / E[M{p,t)]dt 
' j>l 

+ j^^E l-exp(^-^^^s^l|,^>i}M(i'^)(p,5,-t)^ I{ds)dt. 

According to Lemma ll(ii), the first component of this sum is finite and, 
for all r/ > 0, there exists some i.i.d. r.v. D(^^p^ having finite moments of 
all positive orders and independent of (s(tj),i > 1) such that the second 
component is bounded from above by 



Ji 



E 



1-exp -^s 



p-a{:p-l)/{a+r))- 



^-{p-l)/(a+r,) 



/((is) dt. 



Using the fact that /^^ S[l-exp(-t-'^X)] = /o°°(l-exp(-t-'^)) 
for nonnegative r.v. X, one sees that this double integral is equal to 



(i-exp(-r(p-i)/("+'')))(it 



{pri+a) / {a+ri) . 



{a+ri)/(p-l) 



I{ds)E[D 



(1) (a+r,)/ip-l)i 

{v,p) ^ 



If p > 1 + Q + r/, the first integral in this latter product is finite. So, take rj 
small enough so that p > 1 + a + rj and notice then that 



(16) 



{pri+a)/{a+r!)- 



(a+»7)/(p-l) 



< 



{sj>l} 
{pr]+a)/{p-l) 



I{ds) 



^{s,>i}I{d^) 
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The integral (15) is therefore finite as soon as the integral in the right-hand 
side of (16) is finite for some r/ > small enough. Hence, we get (i). 

The same argument shows that Ugtat G for p sufficiently large when 
there exists some e > such that J^i sf l{sj>i}.I(cis) < oo. Indeed, let p > 
1 + a + rj. It suffices then to show that the integral on the left-hand of (16) 
is finite and to do so, we replace the upper bound there by 

„ / \ {"+'?)/(p-i) 

which, by Holder's inequality, is finite as soon as p is large enough and rj 
small enough. 

(ii) We now turn to the proof of assertion (ii) and that Ustat ^ l^'^'^ when 
(H3) holds. The integral (15) is bounded from below by 

roc r 

> / s^" / (l-exp(-4rt-(P-^)/°))P(M(p,t)>rt-(P-^)/°)(it/((is). 
Ji^ Jo 

According to Corollary 3 in [8], the hypothesis (H3) ensures thatt(P-i)/"M(p,f) 
converges in probability to some finite deterministic constant as t — > oo. 
Hence, taking r > small enough and then to large enough, one has P{M{p, t) > 
^■^-{p-i)/") > 1/2 for t>to and, therefore, the integral (15) is bounded from 
below by 



which is infinite as soon as p < 1 + a or J^i s^^^^ ^^l|<,j>jg}I((is) =00. □ 

Proof of Proposition 10. For the standard fragmentation F, let 
^{e,oo)i't) '■= Z]fc>i l{Ffc{t)>£} denote the number of terms larger than e present 
at time t. Under the hypotheses (H3), (H4) and a > — 1, Theorems 4(i) and 
7 of [19] describe the behavior of -/V(£ ^^^(t) as e — > 0. Theorem 4(i) states the 
existence of a random function L such that J2k>i ^kit) = Jt° -^(^) du a.s. for 
all t. Then Theorem 7 says that 

(17) e^+"iV(e,oo) (i) KL{t) as e ^ 

a.s. for almost every t, where K = {I + a) /a^ E[i{l)]. Note that the sum 
Ustat (e) can be written as 

(18) U.tat(e) = E 4A,{t.),oo)(^i 
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are i.i.d. copies of iV(. independent of ((s(tj), tj), i > 

(i) Let T^^'^^ be the first time at which F^*'-'-' reaches 0, i,j >1. With the 
same argmnents as in the proof of Theorem 7(i), one sees that, with probabil- 
ity one, there is at most a finite number of ti < sup j^i{T^^'^^ sj"' {ti)) if and 
only if E[si\pj^iT^^'^^sJ°']I{ds) < oo. This integral is finite by assump- 
tion. A moment of thought then shows that there is at most a finite number 
of integers i,j > 1 — independent of e — such that ^j^l'/g .(^^.^ oo)(^j 
Consequently, the sum (18) involves a finite number of nonzero terms and 

where the functions L^^'-'^'s are i.i.d. and distributed as L. This limit, which 
we denote by X, is null as soon as Ustat = 0, that is, as soon as there is no 
integer z > 1 such that ti < sup j^i{T^^'^^ sj"' (ti)) . This occurs, according to 
the Poissonian construction, with a positive probability. On the other hand, 
the Lebesgue measure of Bl ■= {x > 0:L{x) > 0} [denoted by Leh{BL)] is 
a.s. nonzero and then P{X > 0) > 0. 

(ii) Suppose /^i sj""l|s^<i}/((is) = oo and let Bj^(i,j) := {x > 0:L^^'^\x) > 
0}, which are i.i.d. copies oIBl- One checks that there a.s. exists a time ti G 
Uj>i s]'^{ti)BL(r,3) if and only if the integral /;i S[Leb(Uj>i s7°'Sl(i,j) 

is infinite and that this integral is indeed infinite here, according to the 
assumption on / and since Leb(;BL) > a.s. From this, we deduce that 

^ L^''^\s'^{ti)ti)sY'^{ti) > a.s. for N large enough 

and then, by (17) and (18), that liminfe_o e^+°Ustat(e) > 0. □ 

3. Rate of convergence to the stationary distribution. We are interested 
in the convergence in law to the stationary regime Ugtat- It is already known, 
according to Lemma 5, that, for every random u G P, the process 
converges in law as t ^ oo to the stationary state Ustat, provided it belongs 
to T> a.s. The aim of this section is to strengthen this result by providing 
upper bounds for the rate at which this convergence takes place. The norm 
considered on the set of signed finite measures on D is 

||/i|| := sup 

/1-Lipschitz,supggi, |/(s)|<l 

By / is 1-Lipschitz, we mean that |/(s) — /(s')| < d{s,s') for all s,s' gD. It 
is well known that this norm induces the topology of weak convergence. 

The main results are stated in the following Theorem 12. In case a < 
0, the rate of convergence depends on / and it is worthwhile making the 



Jv 
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result a little more explicit. This is done, under some regular variation type 
hypotheses on /, in Corollary 13. 

Theorem 12. The initial states u considered here are all determinis- 
tic. 

(i) Suppose that a <0 and J^i J2j>i sj°'^{sj>i}^ids) < Then, for ev- 
ery 7 G [l,r] [r is defined by (9)], there exists a positive finite constant A 
such that, for every u satisfying J2j>i exp(— ) < oo, 

||£(F/H(t))-£(U,tat)|| 

= O l^i-(^-i) J2 *7"^ eM-At"'s]^)I{ds) + exp(-Ai^<^) j 

as t ^ oo. 

(ii) Suppose that a = and J^i Y^j>i s^~^'^I{ds) < oo for some e> 0. Then 
for every u E l^~^^ and a < cl){e) /(2 + e) , 

||£(F/(")(t))-£(Ustat)|| =o(exp(-at)) ast^oo. 

(iii) Suppose that a> and /^i J2j>i s^I{ds) < oo for some p> 0. Then, 
for every u£P and every a < 1/a, 

||£(F/(")(t)) - £(U,tat)|| = o(t-°) as t ^ oo. 

Note first that, by Theorems 7, 8 and 9, the assumptions we make on 
/ imply in each case that Ustat G a.s. In case a < 0, the given upper 
bound may be infinite for some 7's. The point is then to find the 7's in 
[l,r] that give the best rate of convergence. This is possible, for example, 
when Jii J2j>i '^{sj>x}I{ds) behaves regularly as x — > 00. In such case the 
statement {i) turns to: 

Corollary 13. Suppose a <0 and fix u such that J2j>i^^P{~'^'j) < 
00. 

(i) // Jii J2j>i ^{sj>x}Iids) ~ l{x)x~^ as X ^ 00 for some slowly varying 
function I and some Q> 0, then, provided —a<Q, 

||£(F/(")(t)) - £(Ustat)|| = 0(/(ti/l"l)t-(^/l"l-i)) as t ^ 00. 

(ii) // — ln{JiiJ2j>i^{sj>x}I{ds)) ~ l{x)x^ as x ^ 00 for some slowly 
varying function I and some q> 0, then there exists a slowly varying func- 
tion I' (which is constant when I is constant) such that 

||£(F/(")(t)) - £(Ustat)|| = 0(t-(r-i) exp(-r(t)t^r/(l"|r+^))) as t ^ 00. 
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In the special case when I[si > a) = for some a > 0, 

- £(Ustat)|| = 0{eM-Bf)) 

for some constant B > 0. 

Proof, (i) First, by integrating by parts and then using, for example, 
Proposition 1.5.10 of [10], one obtains that, for 7 E [1, Q/{—a)), 

I E s7""l{.>.-}/(^s) p. as X - 

(the notation ~ means that the functions are equivalent up to a multiplica- 
tive constant). Then, using Karamata's Abelian-Tauberian theorem (Theo- 
rem 1.7.1' of [10]), one deduces that 

/ y. ^T"' exp(-ts°^)/(ds) PS /(i-i/"7)ii+e/«7 as t ^ cx). 

//I ^ ^ ^ 

Now if —a < Q, statement (i) of Theorem 12 applies and one can plug the 
above equivalence into the upper bound obtained there, hence, the conclu- 
sion. 

(ii) Let 1 < 7 < r. By integrating by parts and then by using Theo- 
rem 4.12.10 in [10], one sees that — ln(/;i '}Zj>i'^^"'^^{sj>x}I{d^)) ^l{x)x^ as 
X — > 00. According to de Bruijn's Abelian-Tauberian theorem 4.12.9 in [10], 
this implies that 

(19) - In^^ ^^s7"^exp(-ts"'^)/(ds)^ ^ f{t) as t ^ 00, 

where f{t) = l/^^{t) with ^(t) = ^{t)/t and $^(t) = f«?/("T)//(iV(-"T)). 
Here *&^(i) = sup{n > 0:(/)(n) > t} and similarly for Therefore, f{t) ~ 
/(t)t^/(^"'"l°l''') for some slowly varying function I (to invert regularly vary- 
ing functions, we refer to Chapter 1.5.7 of [10]) which is constant when 
I is constant. The assumption we have on / allows us to apply Theorem 
12(i) and the conclusion then follows by taking there 7 = T and using the 
equivalence (19). The special case when I(si > a) = is obvious. □ 

Hence, our bounds for the rate of convergence depend significantly on / 
when a < 0, whereas they are essentially independent of / when a > 0. Also, 
in any case they are essentially independent of the initial state u. 

We now turn to the proof of Theorem 12, which relies on a coupling 
method that holds for P-valued X-processes with immigration, as defined 
in Section 2.1. We first explain the method in this general context and then 
make precise calculations for fragmentation with immigration processes. In 
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this latter case, if c, v and / are fixed so that /(si > 1) = and if a varies, 
one sees (without any calculations, just using that particles with mass < 1 
split faster when a is smaller) that the employed method provides a better 
rate of convergence when a is smaller. When /(si > 1) > 0, the comparison 
of rates of convergence as a varies is no longer possible because particles 
with mass larger than 1 split more slowly when a is smaller. 

Proof of Theorem 12. Let X be a P- valued Markov process with 
the superposition property and / an immigration measure such that the pro- 
cesses Xl^^^ ^ u G 2?, defined by formula (11), are P-valued X-processes with 
immigration. Let then ((s(tj), tj), z > 1) be the atoms of a Poisson measure 
with intensity /(ds) dt, t > 0, and suppose that the stationary sum Ugtat 
constructed from ((s(tj), ij), i > 1), as explained in (12), belongs a.s. to T). 
Suppose, moreover, that X(")(t) for ah u e P. 

Then, fix u G P and consider X^") and some versions of X 

starting, respectively, from u and Ugtat- Consider next an X-process 

with immigration starting from 0, independent of X^") and X^^^'^^^ . Then, 
the processes and X/(Ustat)^ defined, respectively, by X/(")(t) := 

X(")(t) +X/W(i) and X/(UBtat)(i) :=X(U='-')(t) +X/W(t), t > 0, are X- 
processes with immigration starting, respectively, from u and Ustat- 

Let now r be a deterministic function and call ri"'' the first time t at 
which xj"^(s) < r(s) for all s > t and, similarly, Tf"^^^ the first time t at 
which xj^''""^(s) < r(s) for all s>t. Of course, the interesting cases are 
ri"^ < oo and Tr^^^^^ < oo a.s. Such cases exist, take, for example, r = 1. 

Our goal is to evaluate the behavior of the norm ||£(X/(")(t)) — £(Ustat)|| 
as t — > oo. To do so, let / :P ^ M denote a 1-Lipschitz function on V such 
that supgg-p |/(s)| < 1. For all t > 0, we construct a function /^(f) from / and 
r(t) by setting 

^//(0)> whensi<r(t), 
/r(t)is;- \/(si,...,Si(,(t))0,0,...), whensi>r(t), 

where i{r(t)) is the unique integer such that Si(r(t)) > '''{t) and Si(r(t))+i < 
r{t). Clearly, as / is 1-Lipschitz and d(s,s') = supj>]^ \sj — s'j\ for s, s G P, 
|/(s) — /r(t)(s)| < r{t) for every s gV and, therefore, 

|i^[/(X/(")(t))-/(U,tat)]| 

(20) = \E[f{XI^^Ht)) - /(X/(u-t)(i))]| 

< 2r{t) + |E[/,,(,)(X/(")(t)) - /,(,)(X/(Ust.O(i))]|. 

The time ri"^ and the function /^j-j) are defined so that, for times t > Tr"\ 
/r(f)(X/(")(t)) takes only into account the masses of particles that are de- 
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scended from immigrated particles, not from u. Therefore, one has 
and, similarly. 

Combined with (20), this gives 

|i^[/(X/H(t))-/(U,tat)]| 

<2r(t) + 2P(r,(") Vr.^*^*^*) >t) 

since supgg-p |/(s)| < 1. This holds for all 1-Lipschitz functions / such that 
supgg-p |/(s)| < 1 and, therefore, 

(21) ||£(X/(")(t))-£(Ustat)|| <2(r(t) + P(T;(") >t) + P(r(^*-*) >t)). 

The point is thus to find a function r such that the above upper bound gives 
the best possible rate of convergence. 

In the rest of this proof, we replace X by an (a, c, z/)-fragmentation pro- 
cess F, in order to make precise calculations. We recall that F^^\t) 
and that the assumptions of Theorem 12 involving / ensure that Ugtat G 
a.s. for all a G M, so that (21) holds for FI^^\ The choice of the function r 
then differs according as a < 0, a = and q > 0. 

Proof of (i). Here we take r = 0. According to the definitions above, ri"^ 
is the first time at which F^") reaches (it may be a priori infinite) and 
^^stat) ^j^^ g^^^ time at which p(Ustat) reaches 0. As recalled in Section 1.1, 
the first time r at which a 1-mass particle reaches is a.s. finite since a < 0. 
By self-similarity, the first time at which a particle with mass m is reduced 
to is distributed as m~°^T. Hence, by definitions of p(") and p(Ustat)^ 

r/")= sup n-^T^j) and r^*^*^*) = sup (s-"(ti)T^*'^) - ti)+, 

where {T^^\j > 1) and {T^^'^\i,j > 1) denote families of i.i.d. copies of r 
such that {T^'^'^\i,j > 1) is independent of ((s(tj), t,), i > 1). 
Now fix 7 G [1, r]. On the one hand, one has 

> tu"), 

J>1 
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which by (8) is bounded from above by X]j>i 6xp(— C;^nn"'^) for some 
constants C^,C;^ > 0. Let < e < C;^. It is easy that this sum is, in turn, 
bounded for all t > 1 by B exp{—{C'^ — e)f'u^^), where is a constant (de- 
pending on 7, £ and u, not on t > 1) which is finite as soon as X]j>i exp(— n") < 
oo. On the other hand, 

P(t^(^*^') >t)< r I Y,P{T>{t + v)s'^)I{ds) dv, 
which, again by (8), is bounded from above by 

for t > 0. Hence, the result. 



We thus have to choose some function r / 0. By assumption, /^i J2j>i s^j'^^I{ds) < 
oo for some e > 0. So, fix such e, fix > 1 and set a := + 7^(1 + e)). 



Proof of (ii). When a = 0, the fragmentation does not reach in general, 
'^e thus have to choose some functii 
) for some e > 0. So, fix such e, f 
Then take r(t) := exp(— at), t>0. 

In order to bound from above P(rr"^ > t) and P(rr^*^*^ > t), introduce, 
for all X > 0, 

Ta,x = sup{t >0:Fi{t) > xexp(-at)} 

the last time t at which the largest fragment of a standard fragmenta- 
tion process F starting from (1,0,...) has a mass larger than a;exp(— at). 
Here we use the convention sup(0) = 0. This time Ta^x is a.s. finite because 
exp(at)Fi(t) when < a < supp>Q as explained in [8]. More pre- 
cisely, one can show the existence of a positive constant C(a) such that 

(22) P{Ta,x >t)< C(a)x"(^+^) exp(-at) for ah x > 0, t > 1. 

Indeed, let t > 1 and note that 

P{rit>Ta,x >t)< P{3ue [t,i]t[:Fi{u)e-xp{au) > x) 

< P{Fi (t) exp(ar?t) > x) (as Fi \ ) 

< x-(i+^) exp(ar/(l + e)t)^[(Fi(t))^+"]. 

This last expectation is bounded from above by -E'Efc>i(-^fc(0)^^^] = 6xp(— (?l)(e)t), 
which yields P{r]t > Ta,x > t) < x~^^^^^ exp(— at), since a = (j){e) — arj{l + e). 
Then, setting C(a) := X]n>i 6xp(— a(r/"'~^ — 1)), one obtains (22). 

By definition, r^"^ is the supremum of times t such that F^^'^ (t) > exp(— at). 
Hence, there exist some independent random variables T^yy^., J ^ 1) where 
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a i/n same distribution as Ta^i/uj^ such that 

T-(u) _ (i) 

Then, by (22), 

(23) P(r,(") >t)< C{a) exp(-at) ^ 

Next, by definition of Tr^^^^\ there exists a family of r.v. Taclp(at )/s (t )^ 
i, J > 1, such that 

^(stat)_ / (ij) + 

and, conditionahy on {{s{ti),U),i > 1), r^^^^uy^^^t^) T„,cxp(atO/«.(t.)' ^'-^^ ^ 
(i i) 

1, and the ''"q exp(at )/s (t independent. This imphes that 

p(4--)>t)<EE^ri:£,(a,)/..(t.)>*^+^) 

i>lj>l 

and then, by (22), that 

C{a 

p /n. 

i>i 

Combining this last inequality with (21) and (23), one obtains 
||£(F/(")(t))-£(Ustat)|| 

< 2exp(-at) [ 1 + C(a) ^ u]+' + (2a)-iC7(a) / ^ s)+^/(ds 
V i>i j>i 

This holds for every r}> 1 and, therefore, 

||£(F/H(i))_£(Ustat)||=0(exp(-at)) 



p( (Stat) > ^) < exp(-at) / y s]+'I{ds). 

2a + e Jii"^^ 



for every a < (p{e)/{2 + e), provided u G Z"'^"''^. Then, as this holds for all values 
of a in an open interval, one can replace 0(exp(— at)) by o(exp(— at)). 

Proof of (iii). Fix < a < 1/a and set r(t) := t~", t > 0. By assumption, 
there exists some p> such that Z]j>i s^I{ds) < oo and we call z the real 
number such that za'^{a + 1) = p{l — aa — az). Note that < z < — a. 
Define then, for x > 0, 

Ta,^:=sup{t>0:Fi(t)>xt-°}. 
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The fact that z G (0,a"^) allows us to choose some 77 > and g > 1 such 
that — aq = q{a~^ — a — z), which, by definition of z, is also equal to 
qza{a -\- l)/p- According to Lemma ll(ii), there exists an r.v. D(j^^q^ with 
positive moments of all orders such that 

a.s. for every t > 0. This implies that 

P{Ta,x >t)<P{3u>t: u'"'F^{u) > x«) 

<P{3u>t: Z)(^,g)n-'?^"("+^)/f > x") 

where B := ElDf^^"'"^] < 00. 

A moment of thought shows that the times r^"^ = sup{t > 0:F^^\t) > 
t-"} and ri'*''*^ = sup{t > : F/"=*"*^(t) > t"'^} satisfy 

r(")=sup(n-"r(^L.i) and ri^'^*) < sup (s-"rHLi - 

where the r.v. r^''\„_i, 7 > 1, are independent such that T^''\a-i , aa~i 
and, conditionally on {{s{ti),ti),i > 1), the r.v. T^^'''}a-i, hj ^ are inde- 
pendent such that r^*'i^a_i ^= s<^<^-^ ■ Using then the upper bound P{Ta,x > 
t) < one obtains 



(a+l) -a(a+l)+p{l-"a)/2o 



which is equal to Bt 5I]j>i by definition of z. Similarly, one obtains 
P{4''''''> >t)< a-^Br" I sp{ds). 

Hence, by (21), 

||£(F/(")(i))-£(U,tat)||<i?t-"fl + E^'+ / T.'WA^ 

V j>i -^'S^i / 

where i? is a finite real number depending on the parameters of the fragmen- 
tation and on a, but not on t and /. This holds for all a G (0, 1/a), which 
gives the bounds o(t"'^), a < 1/ct, claimed in the statement. □ 
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4. An example constructed from a Brownian motion with positive drift. 

Let B he a standard linear Brownian motion and for every d> 0, consider 
the Brownian motion with drift d, 

B(^a){x) ■= B{x) + dx, x>0. 

For any t > 0, define 

L(d) (t) := inf{a; > : B^^^ (x) = t} R^a) (t) := sup{x > : {x) = t} , 

the first and the last hitting times of t by Clearly, < L(^fi-j{t) < 

< ^-s-i since d > 0. It is thus possible to consider the decreas- 
ing rearrangement of lengths of the connected components of 

£(^d){t) ■.= {xe md){t),R(d){t)]:B^a){x) > t}, 
which we denote by F/(^)(t). 

Proposition 14. (i) The process > 0) is a fragmentation 

immigration process with the following parameters: 

• as = -1/2, 

• CB = 0, 

• i^Bisi + S2 < 1) = and z/b(si G dx) = \phr^x~'^l'^(\ - xY'^l'^ dx, x £ 
[1/2,1), 

• I{d){s2 > 0) = and /(rf)(si G dx) = y^{2-ir)-'^x ^^'^ exp{-xd'^ /2) dx, x > 0. 

(ii) The process is stationary. The stationary law is that of a Cox mea- 
sure (that is, a Poisson measure with random intensity) with intensity T{d)\/ (Svr)" 
exp(— dx, x > 0, where T(d) is an exponential r.v. with parameter d. 

(iii) There exists a constant L £ (0,oo) such that, for every u£V satisfy- 

— 1/2 

ing J2j>i^w{~Uj )<oo, an [as, cb^i^b, 1(d)) -fi"'^9''^^''^t(''t'i'0n' immigration 
F/^"^ starting from u converges in law to the stationary distribution £(Ustat) 
at rate 

-£(Ustat)|| =0(t-^exp(-Lt)). 

Note that the immigrating particles arrive one-by-one. 

The fragmentation part of this process, that does not depend on d, is 
a well-known fragmentation process that was first constructed by Bertoin 
in [7]. Let F^'^ denote this fragmentation starting from 1= (/,0, ...). It is 
a binary fragmentation, that is, each particle splits exactly into two pieces, 
which is constructed from a Brownian excursion e^^ conditioned to have 
length / as follows: 

(24) F^\t) := {lengths of connected components of {x G [0, Z] :e^j^{x) > t}}^ 
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for all t > 0. In [7] it is proved that this process is indeed a fragmentation 
process with index = — 1/2, no erosion and a dislocation measure ub as 
given above. 

Proof of Proposition 14. (i) According to Corollaries 1 and 2 in 
[25], the process defined by 

y(rf)(x):=i?(rf)(x + i?(rf)(0)), x>0, 

is a BES'^(3,(i) (which means that it is identical in law to the norm of 
a three-dimensional Brownian motion with drift d) and is independent of 
(i?((i)(x),0 < X < i?(^)(0)). This last process codes the fragmentation of par- 
ticles present at time 0, whereas the process Yf^^-^ codes the immigration and 
fragmentation of immigrated particles. More precisely: 

• Let e^^\ . . . , \ . . . denote the finite excursions of above 0, with 
respective lengths /i,/2, • • • • The Cameron-Martin-Girsanov theorem im- 
plies that the > 1) are the finite jumps of a subordinator with Levy 
measure {87r)~-^x~^^'^e~^'^'^^^ dx, killed at an exponential time with pa- 
rameter d, and that conditionally on > 1) the excursions e^^\e"^ \ . . . 

are independent Brownian excursions with respective lengths li, . . . ,li, 

This gives the distribution of F/(^)(0) = (/i, Z2, • • • ^^'^ implies that the 

process (F/;"^p'^^^'^\t),t>0) defined by 

■^^(d)^''*'^'^^' {lengths of connected 

comp. of {xG [L(rf)(t),i?(rf)(0)]:S(rf)(2;) >*}}■•• 

is an (—1/2, 0, i^b) -fragmentation starting from F/(^)(0). 

• Let J(Y^^^){x) := mfy>xY{d){y)^ x > 0, be the future infimum of y(d)- Oiie 
has to see «^(Y{^) ) as the process coding the arrival of immigrating particles 
and Y'l'^) — -^(Yf^)) as the process coding their fragmentation. According to 
a generalization of Pitman's theorem (Corollary 1 in [25]), ( J(Y(^j), ^(d) — 
•^(^(d))) distributed as (M(d), M(rf) -B(rf)), where M(rf)(x) := sup^^^j %)(y), 
X > 0. Moreover, according to the Cameron-Martin-Girsanov theorem, 
M(^) is distributed as the inverse of a subordinator with Levy measure 

I{d) {si G dx) = ^ {2tt)-^x~^/'^ exp(-X(i^/2) dx, x > 0, 

and, conditionally on their lengths, the excursions above of — 
are Brownian excursions. Let ((A(^)(tj),tj),i > 1) denote the family of 
jump sizes and times of the subordinator inverse of M(^). The sequence 



{lengths of connected 



comp. of {xG [%)(0),%)(t)]:%)(x) >t}} 
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is the decreasing rearrangement of masses of particles that have immi- 
grated at time ti<t with mass A(^)(tj) and that have split independently 
(conditionally on their masses) until time t — ti according to the fragmen- 
tation (— 1/2, 0, z/b). 
• F/(^)(t) is the concatenation of (t) and FI^^!j'''^^^^'°°\t), which 

leads to the result. Note that /(^j satisfies the hypothesis (HI). 

(ii) That F/(^)(t) ^= F/(^)(0) is a simple consequence of the strong 
Markov property of B applied at time Lf^j^^{t). The stationary distribution 
C{FI ^^-^(0)) is calculated in the first part of this proof. 

(iii) It is easy to check that the i/^-dependent parameter Tb [defined 
in (9)] is here equal to 2 and that 

-ln(y I(d)isi & dy) ] — asx^oo. 
Then we conclude with Corollary 13 (ii). □ 

Remark. Let Y(^d) be a BES°(3,(i), d>0, and set 

Fly^j^^it) := {lengths of connected 

comp. of {x E [LY^^^{t),RY^^^{t)]:Y^a)ix) > t}}^ , 

where Ly^^^^ (t) := inf{x > : Y(^) (x) = t} and Ry^a) (0 •= sup{x > : Yj-^) (x) = 
t}. According to the proof above, FIy^^-) is an (— 1/2, 0, z^s, /(rf))-fragmentation 
with immigration starting from (clearly, this is also valid for d = 0). Re- 
call then the construction of the stationary state Ustati as explained in (12). 
It is easy to see that Ustat has the same law as the point measure whose 
atoms are the lengths of the excursions below of the process obtained 
by reflecting Yj-^) at the level of its future infimum. By Corollary 1 in [25], 
this reflected process is a Brownian motion with drift d. Therefore, if d > 0, 
Ustat G a.s. and the stationary distribution is that of the reordering of the 
lengths of the excursions below of a Brownian motion with drift d, which 
is indeed the distribution of F/(^)(0) (by Girsanov's theorem). On the other 
hand, if d = 0, Ustat is clearly not in T> a.s. and then there is no stationary 
distribution [which was already known, according to Theorem 7(ii)]. 

At last, we mention that one can construct in a similar way some fragmen- 
tation with immigration processes from height functions coding continuous 
state branching processes with immigration (as introduced by Lambert [22]). 
This is detailed in [20]. 
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5. The fragmentation with immigration equation. The deterministic coun- 
terpart of the fragmentation with immigration process (a, c, i^, /) is the 
following equation, namely, the fragmentation with immigration equation 
(a, c, u, I) 



dt{iiuf)= x° -cx/'(x) + 

Jo \ JVx 



(E) 



■i>i 



X] 



v{ds) \\ii{dx) 



where > 0) is a family of nonnegative Radon measures on (0, oo). The 
measure f^tidx) corresponds to the average number per unit volume of par- 
ticles with mass in the interval {x,x + dx) at time t. The test-functions / 
belong to C^(0, oo), the set of continuously differentiable functions with com- 
pact support in (0,oo). Note that the hypothesis (HI) implies the finiteness 
of the integral JiiYl,j>i f{sj)Hd^) every /gC^(0,oo). In [2], the sta- 
tionary solution to this equation is studied in the special case when a = 1, 
c = 0, ^{si S dx) = 21|^g[i/2,i]} dx and 1^(51 -|- S2 < 1) = 0, I{s2 > 0) = and 
I{si S dx) =i{x)dx for some measurable function i. Here we investigate 
solutions and stationary solutions to (E) in the general case. 

5.1. Solutions to (E). When / = 0, existence and uniqueness of a solu- 
tion to (E) starting from 6i{dx) are established in Theorem 3 in [18]. More 
precisely, the unique solution to the equation starting from 6i{dx) is given, 
for all t > 0, by 



(25) 



{ritJ)-=E 



k>l 



fe Clio, 00) 



where F is a standard fragmentation process (a,c, z^). Now, we generalize 
this to the case when / 7^ 0. In that aim, we recall that some fragmentation 
with immigration processes starting from u gTZ were introduced in (10). 
Recall also that (p is the Laplace exponent given by (3) and that (p = 4'~4>{^)- 

Proposition 15. Let /xq be a nonnegative Radon measure on (0,oo) and 
let u be a Poisson measure with intensity fiQ. Consider then an (a,c, i^, /)- 
fragmentation with immigration {FI^^\t),t>0), as introduced in (10), and 
define a family of nonnegative measures {fJ-t,t ^ 0) by 



(26) 



{l^t,f):=E 



E/(M"^W) 

k>l 



/gC1(0,oo), />0. 



If one of the three following assertions is satisfied: 
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(Al) a > 0, Jii J2j>i Sjl{ds) < oo and x^Q{dx) < oo, 
(A2) a = 0, JiiT,j>iSj'^(,^)l{sj>i}I{ds) < oo and x'^{^^) fio{dx) < 
oo, 

(A3) a < 0, Jii J2j>i s]'^'^'i-{sj>i}Iids) < oo and x^^" fio{dx) < oo, 

then the measures /it, t > 0, are Radon and the family (/Uj , t > 0) is the unique 
solution to the fragmentation with immigration equation (E) starting from 

Of course, FI^^^ is a "usual" P-valued fragmentation with immigration 
process as soon as /Uo[l,oo) < oo. 

Remarks. 1. Notice that, for ah / G (0, oo), / > 0, 



{l^uf)=E 



i>l k>l 



-ti<tj>lk>l 

where {{s(ti),ti),i > 1) [resp. {ui,i > 1)] are the atoms of a Poisson mea- 
sure with intensity I{ds)dt (resp. /^o) and F is an (a, c, i/) -fragmentation, 
independent of these Poisson measures. By (6), this can be written as 



il^tJ) 
(27) 



i?[/(xexp(-C(p(x-t))))exp(e(/9(x-t)))]/io((ix) 



+ / / J2E[f{s,eM-C{p{s'^umeMC{p{s'^u)))]I{ds)du, 



where ^ is a subordinator with Laplace exponent cf). It is not hard to see that 
there exist some dislocation measures ui ^ V2 that lead to the same (f). In 
this case, the previous formula shows that the (a, c, z^i, I)- and {a,c^V2-,I)- 
fragmentation with immigration equations have identical solutions. 

2. Assume that one of the assertions (Al), (A2) and (A3) is satisfied, so 
that the measures pt,t>0, are Radon. Then, these measures are hydrody- 
namic limits of fragmentation with immigration processes. Indeed, let u^") 
be a Poisson measure with intensity n/^o and call F/^") a fragmentation with 
immigration process with parameters {a,c,u,nl) starting from u^"^ Then, 
for every t > 0, 



-F/(")(t) 
n 



vaguely 



Ptidx) 



a.s. 



This holds because F/(")(t) is the sum of n i.i.d. point measures distributed 
as F/^"'^^^(t) for some (a, c, z/, /)-fragmentation with immigration 
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The strong law of large numbers then implies that, for every / G C^(0, oo), 



" k>i 



k>l 



and the conclusion follows by inverting the order of "for every / G C^(0, oo)" 
and "a.s.," which can be done, for example, as in the proof of Corollary 5 
of [18]. 

Proof of Proposition 15. Let /ii,t>0, be defined by (27) [equiva- 
lently, (26)]. 

• It is easily seen that these measures are Radon if (Al) holds. To prove 
this is also valid for assertions (A2) or (A3), we need to evaluate the rate of 
convergence to of P{a < x exp{—(,{p{x°'t))) < 6) as x — > oo, < a < b < oo, 
when a < 0. First, note that this probability is bounded from above by 
P{x ex.p{—(^(p(x"t))) < b), where ^ = '^l{^<oo} is a subordinator with Laplace 
exponent (j) = (p — (/'(O). Then for u>0 and v > 0, 

P{^{u) >v)<{l- e~^y^E[l - eM-v~^'^{u))] 

(28) 11 - 

= (l-e-i)"Hl-exp(-7x0(^;-i))). 

When a = 0, this implies that 

(29) P(a<xexp(-^(t)) <6) = 0(^((ln2;)~^)) as x ^ oo. 

When a < 0, by the definition of p and conditionally on 2x°^t < p{x'^t) < oo, 

2x"texp(a^(2x"t)) < / exp{aCir))dr< exp(a^(r)) dr = 

Jo Jo 

and, consequently, P{2x°^t < p{x"t) < oo) < P(exp(a^(2x"t)) < 1/2) which, 
by (28), is an 0(x") as x — > oo. Moreover, again by (28), P(xexp(— ^(2x"t)) < 
b) = 0{x°) and, therefore, 

(30) P{a < X exp(-^(p(x°t)) <b) = 0(x°) as x ^ oo 
since 

P{a < xexp(-^(/)(x"t)) < b) 

< P(2x"t < p(x"t) < oo) + P(xexp(-^(2x"t)) < b). 

Now, suppose that (A2) or (A3) holds and take /(x) = xl^xG{a,b)}^ < a < 
6 < oo. Using the results (29) and (30), one sees that (/^t,/) is finite. Hence, 
pt is Radon. 
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• Suppose that (Al), (A2) or (A3) holds, so that the measures fit,t > 0, 
are Radon. Consider then the measures rjt, t>0, introduced in (25). One 
checks that 

{fJ't,f)= {r]x'-t,fx)fJ'o{dx)+ y2{7hfu,fs,)I{ds)du, 
Jo Jo Ji^J^i 

where fx'-U^ f{xy), x G (0, oo), / G (0, oo). Theorem 3 in [18] states that 
{ilt^ty 0) is a solution to (E) when / = 0, that is, 

(^t, /) = /(!)+ f\nu,Af)dv, 







where 



(31) Af{x)=x^i-cxf{x) + 



J2f{xs,)-fix) 
■i>i 



This equation rehes on the fact that, for / G C^(0,cxo), y4(id x f){x) = 
x^~^°G{f){x), where G is the infinitesimal generator of the process exp(— ^) 
(see the proof of Theorem 3 in [18] for details). 
Using then that x'^Afx = {Af)x, one obtains 

(32) ivx^u fx) = fix) + tivx'^v, {Af)x) dv 

and, therefore, by Fubini's theorem. 







{l^tJ) = {f^o,f)+ / {i]x'^u,{Af)x)^o{dx)du 



JO 



+ I [r I T.(Vsfv,{Af)s,)I{ds)dv+ I Y.f{s,)I{ds 

(/xo,/)+ f\fiu,Af)du + t f ^/(s,)/(ds). 
Jo Jl^ TC. 



(to see why Fubini's theorem holds, call [a, b] the support of / and suppose 
/ > 0. The same argument holds for the integral involving /). Hence, (/it, t > 
0) is indeed a solution to (E). It remains to prove the uniqueness. This can 
be done with some minor changes by adapting the proof of uniqueness of 
a solution to (E) when 7 = (see the third part of the proof of Theorem 3 
in [18]). □ 

5.2. Stationary solutions to (E). As in the stochastic case, we are inter- 
ested in the existence of a stationary regime. We say that a Radon measure 
/istat is a stationary solution to (E) if the family (//i = //stat ,t>0) is a solu- 
tion to (E). 
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Proposition 16. (i) There is a stationary solution to (E) as soon as 
Jii J2j>i Sj^(.ds) < oo and, conversely, provided that hypothesis (H2) holds, 
there is no stationary solution to (E) when JiiJ2j>i'^j^{ds) =00. In case 
Jii J2j>i Sjl{ds) < 00, the stationary solution figtat is unique and given by 

Mstat(da;) := x""" fi''^°^\dx) , x>0, 

where the measure /^slj'°™^ is independent of a and is constructed from c, v 
and I by 

(33) {f^^:T\f) := r / ^i?[/(.,exp(-^(t)))exp(e(t))]/(ds)di, 

/GCi(0,oo). 

(ii) Suppose /^i J2j>i Sjl{ds) < 00 and xfio{dx) < 00 and let {fit, t'>0) 
be the solution to (E) starting from /iQ. Then, 

vaguely 

Mi Mstat as t^ 00. 

Remarks. 1. It /istat exists, then Ustat G a.s. and (ustatj) = -E^[(Ustat, /)], 
/ G C^(0, 00). Note that it is possible that Ugtat G T^\T^, which then imphes 
that there exists no stationary solution in the stochastic case, although there 
is one in the deterministic case. Conversely, Ugtat niay belong to T> a.s., even 
if its "expectation" measure fi defined by {fi, f) := E[{\Jstat, f)] is not Radon. 
Then there exists a stationary solution in the stochastic case, but not in the 
deterministic one. 

2. Call A := sup{A : /^i J2j>i Sjl{ds) < 00} and suppose A > 1. Then the 

statement (i) and the relations £^[6""^^^*^] = e~*'^^''\ t,q>0, imply that, for 
all 1 + Q < A < A + a, 

(34) / xVstat(fix) = </.(A - a - 1)- M ^ s,^-"/(ds), 

and that this integral is infinite as soon asA>A + aorA<l + a, provided 
0(0) = [which is equivalent to c = z^(X]j>i Sj < 1) = 0]. This characterizes 
//stat and is more explicit than (33). 

As an example, it allows us to obtain the more convenient expression 



00 



Aistat((ix) = (^x '^i{x)+2x " j yi{y)dyjdx 

in case ly is binary, i^(si G dx) = 21{2.g[i/2,i]} dx, c = 0, and /(si G dx) = 
i{x) dx, I{s2 > 0) = {a £ M). This latter result is proved in a different way 
in [2]. 

Others examples are given by the equations corresponding to the fragmen- 
tation with immigration processes constructed from Brownian motions with 
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drift d> (Section 4). The immigration measure /(^^ satisfies /^i J2j>i Sj^{d){d^ 
oo for all A > 1/2 and, therefore, there exists a stationary solution to the 
equation. One can use (34) to obtain 

fJ'sta.tidx) = — - exp{-x(f/2) dx, x>0. 

This can also be shown by using remark 1 above and the stationary law 
£(Ustat) given in Proposition 14(ii). 

Proof of Proposition 16. (i) We first suppose that there exists a 
stationary solution nt = Ustat, i > 0, to (E). Of course, then dt{fJ.t, /) = for 
every t>0 and / G C^(0, oo), and, consequently. 



(/istat,A/) = - / ^/(s,)/(ds). 



where Af is given by (31). Letting t — > cx) in (32), we get by dominated 
convergence that {rjxo't, fx) — > and then that f{x) = — J^{r]x'^v, {Af)x) dv, 
x G (0, oo). Hence, 



(^stat,^/) 



LT. {ris^vMf)s,)dvI{ds). 



We point out that this formula characterizes //gtatj since A(id x f){x) = 
x^~^'^G{f){x), where G is the infinitesimal generator of exp(— ^) and since 
G(C^(0, oo)) is dense in the set of continuous functions on (0, oo) that vanish 
at and oo. Using then the definition of rjt and formula (6), one sees that, 
for every measurable function g with compact support in (0,oo), 



(35) 



p POO 

/ E^7"/ E[9isjexp{-av)))eMa + c^nv))]dvIids), 

i>i 



using for the last equality the change of variables v p{s"v) and that 
exp(a,^p(t,)) dp(u) = dv on [0,D), D = mi{v: (,p(^y) = oo}. This gives the re- 
quired expression for fistat- 

Note now that the previous argument implies that a stationary solution 
exists if and only if 

E[g{sj exp{—(^{v))) exp{£^{v))] dvl{ds) < oo 
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for all functions g of type g{x) = xl^^^^^^i,^, < a <b. For such function g, 
the previous integral is equal to 

(36) / E ^A{s.>a}E[Ti^s,/a) - Tt.^s,/,)]Iids), 

where t/ := inf{u:^(ii) >t},t> 0. If hypothesis (H2) holds and ^ is arith- 
metic [i.e., if (H3) holds], the renewal theorem applies (see, e.g.. Theo- 
rem 1.21 in [5]) and -^[Tj^^^y^^ — converges as t — > oo to some fi- 
nite nonzero limit. In such case, the integral (36) is finite if and only if 
Jli J2j>i Sj^{sj>i}I{ds) < oo, y b > a > 0, and, therefore, there exists a sta- 
tionary solution if and only if JiiJ2j>iSj'i-{sj>i}^{ds) < oo. This conclu- 
sion remains valid if (112) holds and ^ is not arithmetic, since the renewal 
theory then implies that limsup^^g^ E'fTj^^^^^^ ~'^in+{t/b)^ < oo, and that 
liminft_^oo ~ ^^L+it/b)^ > as soon as ln6 — Ina is large enough. 

Last, to conclude when (112) does not hold, remark first that = T^ Ae{k) 
[the subordinator ^ and the exponential r.v. e(A;) are those defined in Sec- 
tion 1.1] and then that 

ElT^f , . - r5 + , „J < E\t}, , , - t}., < E\t} J < oo. 

L \n(sj/a) ln+(sj/fe)J — L \n(Sj/a) ln+(sj/6)J — ^ ln(6/a)J 

In this case, the integral (36) is finite as soon as /^i X]j>i ^j'^{sj>i}I{d^) < c«) 
V6>a>0. 

(ii) Under the assumptions of the statement, the measures ^t-, i > 0, are 
Radon and therefore satisfy (27) for all continuous function / with compact 
support in (0, oo). The integral involving /xq converges to as i ^ oo, since, 
with the assumption xfj,o{dx) < oo, the dominated convergence theorem 
applies. Hence, (/it,/) (/^stat,/), using the definition (35) of //stat- □ 

t — ^oo 
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